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i= “OBWESTIVES 


The objectives of the project were: (1) to continue the 
visual, or graphical, analysis o£ the BRIGHAY data (the 
Miberal GitOoLe of Which was “reportad in Ref. 1) in an 
attempt to verify anomalous signal behavidr aiready 
documented and to discovec any new anomalies whica might be 
present; (2) toward that end, to make whatever changes were 
necessary in both the programs and data nandling procedures 
to meet either new systen constraints or to enhante over-all 
pEQ}Ject efficiency; and (3} to analyze the data from a 
statistical Standpoint (via two entirely new computer 
programs) in an etfort ts dateraine how often ani by what 
amount (in decibels) s2veral different classes of signals 


change. 


Ene sw LOLWdt ~O@L th= data, particularly that of the 2.3 
Minute time duration limitation (although theres w3re several 
time-wise contiguous data sets), nade long t2ra analysis, 
both visual and statistical, impossible. In addition, the 
dally intercept times for the data sets precluded any 
Bxgnificant effort b2iag made to correlate Slgnal 
characteristics from even approximately the same time of the 


day. 





Ei INTRODUCTION 


It is well documentei that as radio signals tcaverse the 
ionosphere (including signals bounced off of on2@ of its 
layers) they are often cnanged in a manner that joes beyond 
Simple attenuation due t> the distance travelled. Both the 
make-up £ the ionosphere itself (how many ions are present 
at any one tine, the elusive Sporadic-& O9nenomenon, 
night-day asymmetry and the geomagnetic non-reciprocity to 
name but a few) and the more, dramatic disruptions of 
lonospheric behavior (solar flare caused magnetic storn 
activity, polar cap absorption and th2 auroral displays are 
three) team with the aore common, every day paenomena of 
"“Paraday Rotation" and asisea-producing lightning storms co 
alter any signal with which any or all of ta2m come in 
contact. The results caa p2 severe attenuation, pronounced 
amplification, "ducting", "whistlers", and fading (Sometines 
Satled the Dellinger srfect). References 2 - 6& were 


excellent sources for ionospheric phenomena infornation. 


PeoWosomecrtewbat the application of graphical technicguas 
to signals waoich had traversed the ionosphere aigat add some 
Mmosignt into the problem of exactly what happens to these 
Signals. 


For the most part, th2 visual prasentation techniques, 
ShoWn in various stages in Fig. 1, used in this study ware 
the same as thos2 presanted in kef. 1. HoWsTer, sone 
Changes have been made and will be discussed in detail in 
the appropriate section. As in Ref. 1, interestiag sets of 
scans have been reproduced as graph plots (using the 


mamoatThC plotter) for Eurther analysis and presentation in 
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this paper. All thirty-four of the previously wunanalyzedi 
complete data sets wer2 studied with varying degrees of 
success concerning anomalies (or, least, percaeivei to be sod) 


found therin. 


The swept-tuned receiver and associated eguipment used 
to acquire and record the data sets is an Operational 
receiver eguipment used for other than experinental data 
gathering. Its proper operation regaires frequent tuning of 
numerous senSitivity adjustments. How often ani by what 
amount these attenuators and thresholds are changed have a 
profound effect on the performance of the equipment. 
Therefore, programs MIN/4AX 1 and MIN/MAX 2 were written in 
meer roOort tO accumulate 215ugh Statistical data from the 
Sets to warrant the rcecommeniation of some appropriate 
pegGoOrEIthm@ —£Or  Controlliag the attenuator aljustments. 


sectlom ee LY Of Chis paper will addrass itself td that topic. 


A. NATURE OF THE DATA 


The data was gathered using a Sw2pt-tuned rec2ziver with 
either an omnidirectional or 12 degre= beam antenna. The 
data sets were produced oy sampling the receiver output at a 
25 HZ sampling rate as it scanned Towieato. CULOUsa an ooo 
KHZ band. Each data s2t ztonsists of, at most, 3597 scans 
Meach scan being compss2a of 318 contiguous 2.8 kHz wide 
bins moving down in freguancy as the bin numb2r increases) 
Peaced On /194 sards. Caici one gives anplitude laformation, 
Memeecrarbenign Guahtization levels, for bins 1 through 150 
Peele 10LOLMatioOn On Dins 161 through 318 is found on card 
two. The quantization levels are given in 5-odit binary 
words with even numbered bin information found in rows 11 
emmougr 3 while oid numb2aced bin infocmation is in rows 5 


through 9. Pema cddikerog, ach Card has the card aumber ani 
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set number punched in row 12 for accounting purposes. With 
each set there is a paranster card which delineates the date 
and time of Signal acgaiisition, base frequency, antenna 


configuration and attenuation applied. 


Sumppecoxumately lo) data sets cecorded, forty-four were 
made available for this study. Figure 2 is a table of 
pertinent information conserning the forty-two conplete sats 


analyzed in Ref. 1 and this paper. 


Be CHANGES IN THE GRAPHICAL ANALYSIS PROCEDURE 


Equipment reconfigurations sins2 the graphical analysis 
procedure was formulat2i and reported in Rer. 1 and efforts 
to make the data handling more efficient demnand2d several 
Changes both in proceliur2 and in th2 two nain analysis 
programs, Anomaly A and Anomaly 8. Sosecifis Latormacion 
Soncernaing the workings 9f the two programs may b2 found in 


Pam Darra culariy in section II 8. 


First of all, the sheer number Of binary cards to be 
transcribed onto magnetic tape (nearly a quarter of a 
million) obviated the card-to-tape Pos) Ol que used 
previously. A much faster and far noce celiaole procedure 
WwaS initiated wher2in the Bi 360/67 system, belonging to 
the W. R. Church Computer Center, was used after nornal 
WOrking hours. The cards were first cead onto 9-track tapes 
(800 BPI, even parity) and then transcribed to the J7-track 
tapes (556 BPI, odd pacity) ESewered LOE the 795-9300 
computer. Partial dumps were made of all sets from both the 
9- and 7-track tapes to ensure that the cards r2ai ontod the 
Capes properly ani that ta2 tape-to-tape translation w#as 
Successful. The programs used to read the cards onto tapes, 


translate the 9- to 7-track and dump the tapes foc 
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G@uicilttyeGontLol Can be £5und at the end 3£ this paper. 


This new system of data loading resultiag in the 
requirement to change the sections of the Anomaly A _ and 
Anomaly B programs refering to data read-in. This was 


HecOomplished with little dierficulty. 


Next, a metasymbol subprogram was added to 2ach of thea 
analysis programs to enable entire files (data sets) to be 
skipped when a spacific file requir2i was locat21 other than 
eirst Of a tape. Beawscdiiang | the  coucin=  "FORSCN'’ 4 
kKeyboard-entered number of files (NFILE = ) would be passed 
over. This proved to be axtremaly valuaole as a time and 
labor saving device. Th2 coutin2e "BAKSCN" was alsd available 
MOmmDaAcKINg OVer entire files but cnronis tape drive 
malfunctions negated the use of this even ndcr2e valuable 
tool. 


ene VrINndlestagjOb Os0gram Change was distated by 4 
computer center equipmant resonfiguration wherein the 
CALCOMP 563 plotter was no longer available. Th2 new systen 
tO be used was the PDP-11 computer in conjunction with a 
VERSATEC plotter. [This cnange required modifications of the 
analysis program sections iealing with the processing of the 
Sutput data (desired plots from the AGI-10 screen), Found at 
the very end of the main program, and the generation of a 
new subroutine , cailei the Headec Subroutine, which packs 
the header data into the proper format for the piotting 
Unit. In addition, 2 PDP-11 plot package (used to interface 
the data tape with the plotter) had to b2 extensively 
modified to meet the demands of this anomaly data. Listings 
of the revised Anomaly A and 6B prograns, the aetasymbol 
Peogram for calling "“FORSCN" ani "BAKSCN" and the PDP-11 
Peeteroulcine (Called "aa.c" for Anomaly A plots - this 
coOutine alone is iisted since tha one utilized for the 


mmeomary 6 Graph plots differs only in a faw of the 
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Mooromilateiyerts0 Linesseoapraising “"aa-c") can oe found at 
the back of this paper. Subroutines sSINP, GINPUT , FNS and 
VCD were not changed and, hence, ace not listed in this 


paper but may be found ia Ref. 1. 
oR ANOMALY A AND B PLO?S5 


The plots which have resulted fErom the graphical 
Aamalysis portion of this study reguire some DELOE 
explanation. Due to th2 data handling restrictions imposed 
Eyeetne X0S-9300 Computer, Only a Small portion of the data 
in each set couli be viewed at any one time. This resulted 
mee some dirtficulty in visualizing what could be long-tern 
(time wise) or wide (freqguency wise) phenomena aad demanded 
that each set be analyzed several tines before the entire 


set had been viewed. 


DiesAgoOMaly A (plOES “CONSiSE of twenty scans (signal 
Wavefronts sampled 0.04 seconds apact - or further if scans 
were deliperately skipped to view long-term effacts) with 
mec, mMOSt cvecent scan ia tine falling at the bottom of the 
Pecture. The scans couli be brought onto the scr2en in any 
number up to, and including, twenty at a time. 4s Shown in 
Pegs) 3, the vertical axis iS marked sff in 5 418 increments 
Which apply only to tans bottom scan (called "base scan"). 
The horizontal axis is divided ints ten 28 kHz wide sections 
each composed of ten [rejguency bins 2.8 kHz wide for a total 
ea 100 £requency bins ani 230 kHz.- This width cana be varied 
t> other than 100 bins 1€ desired. The leftmost din (ternei 
the "base bin") is selected by keyborad input, thus allowing 
the analyst HO. sliis up and jown the data set 
(frequencv-wise) at will. The scale along the horizontal 
axis, then, consists of iacscrements of ten bins starting with 


the "base bin" number at the origin. The term "oase_ freg!! 
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} 


refers to the base freyuency of the data set ani decreases 
in 2.8 kHz increments as the bi2 number increases. Shading 
has been provid2d on some of the plots to highlight the 


suggested anomaly. 


The Anomaly B plots are a combination of th2 Anomaly A 
plot technique ani a modification which enables the analyst 
to view selected bins (signals) in the time vs. amplitude 
domain. At tha bottom of each plot are five scans (Anomaly 
A style) used to aid in the selection of the bins to be 
plotted in the upper traces (see Fig. 4). Piers teliron, t5 
the restricted number of scans td be seen at one time (only 
eye, cOmpiuter Limitations reguirei that the aumber 2a 
frequency bins presentei be reduced to eignty. Other than 
those changes, discussio1 3f the Anomaly A plots is also 
valid for the lower set of Anomaly 8 plots. Th2 upper set 
of Anomaly B plots is a sollection of five keybsard-seisected 
Signals chosen from anong the scans below. Phe traces are 
numbered one to five, fron top to bottonu respectively. [he 
MeGcizontal axis 2S a time line composed of aighty data 
Beints at least 9.08 sesonds apart (dependiag upon the 
humber of scnas skippei, if any) with the oldest point at 
the extreme right side of the trace) with the hasa marks on 


the vertical axis refering to 10 1B amplitude levals. 


AS seen to some extent in the Anomaly A plots and, more 
markedly, in the Anomaly B plots, the scans look fracturai 
in appearance. [his was not the cesult of intentional 
programming, ratner alee Was the product of the new 
Pees lii/VeRSATEC plot routine conbination. 
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Figure 3 - SAMPLE ANOMALY A PLOT 
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Pei othPHbGal ANALYSIS 


A. DISCUSSION 


PEVOrW CO saultiatzon of the graphical analysis process 
Via programs Anomaly A ani B, an enumeration of many of the 
possibie "anomalies" expected was made. This was not to 
imply either that Signals which had been markedly affected 
by the ionosphere would only exhibit these graphical 
manifestations of anomalous performance or that ali forms of 
Signal anomalies would be seen by the analyst. Topas 
enumeration procelure was w@arely an attenapt tds ‘'visualize" 
what effects those sigial-changinj ionospheric phenomena 
discussed in the introduction might have on th2 data sets 
yet to be analyzed. Although no attenpt was made to make 
this study's "anomalies" sonform with those ideatified in 
Ref. 1, there has Deen ao attempt to verify any previously 


seen. 


The graphical effects on the signal wavefornas expected 


Mele. 


1. a sudden enhancenenat/depression of one oc more 


Sigudis (tong ob) Shore term) ; 


oo ai DIeeee de ssantinvity ~across an sntire 


wavefront Or portion thereof; 


Seeoeeolgnal Ybugie" or "valley" moving either up 


or down in frequency as the wavefront 


us) 





progresses (Fig. 4, Ref. 1 refers) ; 


(OSG hTEt ta bregqucsncy Of Ona Or more signals 
(either short tern, long term or aperiodic in 


nature) ; 


5. ah unexdlainabl2 enhancement of one or more 
signals with a concomitant depression of one 
or more Signals at the same time in th2 same 


data set; 


6. a complete loss of continuity of tne signal 


wavefront for vaclable periods of tima; 


7. unexpected variations in one or mor2 signal 


GONntours as the wavefront progresses; 


Seman yamoue-oOt-the=- ordinary graphical implication 


of anomalous percformance. 


ESecetcmilr= ert 6 5, 400 Gewere discerned 1n the present 
study while several of these, plus effects 3, 5, and 7, were 


reported in Ref. 1. 


Be. ANOMALY A ANOMALIES 


As shown in Ref. 1 and verified in this paper, most of 
the data sets (when viewed on the graphics terminal) 
exhibited a great measure of uniformity in signal makeup in 
the time domain and ii] not yield any obvious graphical 
maomeliies. in fact, fully fifty per cent of tne BRIGHAM 
data sets failei to provii2 any reportable irregularities in 
Proniat pattern with the other fifty per cent giving fairly 


uniform signal wavefronts for most of the tine recorded. 


20 





The latter group did, however, exhibit anomalous (or, at 
least, assumed to be so) characteristics which were 
verifications of performaas2 reported in Ref. 1 or were 
newly discovered. The aniformity mentioned above aided 
immeasureably in the detection of anodmalies sinzte even a 
Slight wavefront variation could be discerned i1 the midst 


of otherwise uniform signal progression. 


Figures 5 and 5 are two examples of minor 
discontinuities (affectiny 30 t> 40 bins - approximately 100 
kHz) which appear to be the result of a multi-signal 
amplitude enhanc2ment for a very short period of time. This 
effect was founi scattered through almost half of the 
thirty-four data sets vi2wei via Anomaly A. Figuc2 7 is an 
Sxanpke Of a@perigdic discontinuities which ran accoss almost 
ali of the entire data sat bandwidth. This chang=2 did not 
occur over the entire time duration of th2 s2t and, as 
mentioned above, could not be traced over its entire 
frequency spectrun. Pits weuypes Of discontinuity was aiso 
found in data sets 33 ani 80. Both wOrerenhe GLSSOntlouL ey 
examples were anticipated although tae latter forn holds the 


most promise, anodmaly-wise. 


Fading was another effect expected but tie examples 
found were somewhat sSurpcising in that it was anticipated 
that the entire fregueacy bandwiith of a signal wavefront 
(690 KHZ) would be affected. fhis pcovei not co be the 
Bascee wlth dit of the examples of fading implicating at most 
20 frequency bins at a tine. Figur2 8 shows the signal 
Bentered on bin i48 fading in and dat, a process which 
lasted several seconds. Figure 9 shows a much smaller 


time-duration fading centered on bin 24. 
A third typ2 of possiol2 anomalous activity is that of 


mies Shitting, we or down in frequency, of eith2r a single 


Signal or entire 390 kHz wiie signal scans. The orecise 
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Base Scan 141, Base Freq 3.5 MHz, Base Bin l 


Figure 5 - NARROW DISCONTINUITY EXAMPLE #1 
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Figure 6 - NARROW DISCONTINUITY EXAMPLE #2 
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TIME 





FREQUENCY 


set 53, Base Scan 544, Base Freq 6.7 MHz, Base Bin 100 


Figure 8 - FADING EXAMPLE #1 


jE 








FREQUENCY 


Set 33, Base Scan 270, Base Freq 5 MHz, Base Bin 1 


Figure 9 —- FADING EXAMPLE #2 
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explanation of this phenomenon is not available but it must 
be kept in mini that nalfunctions in the recaiver systen 
itself could be the cause rather taan sone 
ionosphere-relatei activity. Figuc2e 10 shows a possible 
iment shitt (approximately 3 bins) and d@1stortion of two 
Side-by-side signals. Pigures 11 and 12 are 2xamples of 
apparent frequency shifts of 160 bins poth of wnich turned 
out to be data set probi2ns. Since each scan is tomposed of 
two data cards of binary coijied amplitudes (bins 1 - 160 on 
card one and bins 161 - 318 on card two), the ori2r in which 
the two cards appear and, as importantly, the existence of 
both cards for 2ach scan are absolutely vital to the correct 
representation of the signal wavefront. For example, if 
card one of 2 Scan iS nissing, than card one oF the next 
scan is taken as card two 9f the former. This has a domino 
effect which will graphically manifest itself as an apparent 
Shift of 160 bins in ail subsequent scans. Figure 11 
MmIiuSsttates precisely tais effect. A manual seacch of the 
missing from the set. Fijiare 12 is an interestinj twist on 
bois problem 12 that the slgnal apparentiy shift2d 160 bins 
Po ethe fight (Or, down ia Ereguency) and then, taree scans 
eter, Shifted rigat back again. Another se@arch of thea 
cards showed that ,as 2xpecsted, two cards were missing, one 
fo Cause the anitial 160 bin fright shift and a second to 
returo the data cards to their correct sejguence three scans 
jater. One last example of signal shifting 1s by £ar the 
most interesting of the lot. As illustrated by Fig. 13, 
there is a shift down in freguency of approxinately 12 bins 
(34 kHz) affecting the entir2 318 bins (890 KHZ bandwidth) 
of aperiodically spacei scans. This phenomenon was s22n 
over the entire length of the data set at apparently randon 
intervals. Data sets 138 and 143 aiso exhibited this 
migeclodic frequency shifting, although the number of bins by 
Which the scans were shifted varied slightly ovec the three 
sets. 
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Set 3, Base Scan 336, Base Freq 3.2 MHz, Base Bin 100 


sc Se Oe oH hPTiNG EXAMPLE #1 


28 








BENG SEL EERIE SPY TB 
ms § Gaaeaereac |» 3 
SESSA Sosa? | 
SS ESE esssissse [ § 


Cae € — —, ~ 
SEPES SECS MOSSES 





Base Bin I] 


Base Freq 18 MHz, 


SHIFTING EXAMPLE #2 


TIME 









jes 
was 
ieee 
/] : 
My v1 
: L 
a 


N 
aN 


AWA 
ANA 

v 

aN 


\ 


V 


Wwe \ 
ee 
AMOS 
AA 
‘ W ye . 
AV 
\ 
AY 
AG 

















= | EC CL RA OR 
SEES ESS Lear OT SS 
SSNS Se 
SESS (SS. 
SS SESE, Ue CEES Yi 
3s <S SF psc Sa55E5. 
Sekikarer Seeerests 
OER CRE 
xs 


os j 
Ss S 
, I SRSEESERS 
— Ee 
cre pal S S58S5S £55 


Cire 
rca Ole Bi eau 





Base Sean es 10> 
Figure ll - 


Set 139, 





TIME 


He we 
ie | i‘ Ae 


NS 


Ah 
orp 
ACOs 


h 
ie 


=“ 
\ 
Ay 


aA 


BRUAVANy 


A ee 


oeeee 
stare eet 


set ate tele’ 


es eet C ERE 


SSS 


30 





Eee a an ee 


FREQUENCY 
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Base Scan 2940, 


Set 126, 
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One fee category for possible signal anomalies 
consists of five differeat phenomena each of whica occured, 
with one exception, only once in the thirty-four data sets. 
Figure 14 is an example wherein three different "anomalies" 
were present. (It should be noted that base scan = 
"various" for Figs. 14 ani 15 means that each 9f the three 
sections of the waveforms in the two figures haji different 
base scans and inaplies that this in no way detracts from the 
usefulness of the figures as examples of the affects in 
gestion.) Bins 1 througa 40 and 83 through 98 exhibit 4 
very high degree of signal attenuation not seen in normal 
data sets (i.e. data sets not found to have either a 
receiver malfunction orf test signal present). Bins 41 
through 87 exhibit normal signal characteristics. Bins 101 
through 200 evidence a combinatida of unusual attenuation 
Mia normal Signal activity, “while bins 201 througa 270 show 
either a general distortion of the signals present or 4 
falrly high amplitude level of anbient noise and/or 
distorted signals plus noise. None of these phenomena could 
be explained by receiver nalfunction, data card sequencing 
or special signal characteristics (1.e. test signal, etc.). 
megure to ShOwS a@ data set exhibiting no csontinuity o£ 
Signal whatsoever. Non-ee Or eeene Signal  “banps’ “can be 
followed in a uniform maaner as time progresses. fhe theory 
that a Serious receiver system malfunction was t212 cause oz 
mire GiStOLtLOR #2aS borne out by a cali to ‘the cognizant 
a@chority. This malfunction was also verified for data sets 
25 and 54. The sawtooth 2ffect seen at the midile, right 
hand side of Fig. 16 is a pattern not seen elsewiere either 
in this set or any of the other thirty-three sets analyzed. 
It shows wild fluctuations in saven adjacent amplitude 
points (2.8 kHz apart) i2 one scan, whereas the crest of the 
scan appears t9 be normal. NOMMexptanation “for this 
paenomenon can be offerei. Another unexpected anomaly was 


the two flat-topped signals, centered at bins 124 and 128, 
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FREQUENCY 


met, o, Base Scan 1815, Base Freq 3.2 MHz, Base Bin 100 


Figure 16 - SAWTOOTH PATTERN EXAMPLE 
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PouUngeeti orig. t7. It appears to be the result of some 
special attenuation appliei to those frequency bins alone, 
because, aS can be seen by the rather "large" signal 
centered on bin 178, other signals appear not to be 
affected. The additional attenuation applied to this set 
during reception , 10 dB, is lower than levels applied to 
other sets which wer? not Similarly affected. The best 
phenomenon observed in tais final category is found in Fig. 
18. For some reason (other than th2 possidility of a 
different form of receivec malfunction than seen defore) an 
Otherwise normal data seat lost total Signal continuity 
across the entire frequency bandwidth with subseqient scans 
peo p1timg 9d higmayedistotted mature. While this could not 
have been caused py a cardi sequencing problem, its source 


Cannot be identified. 
Cs ANOMALY B ANOMALIES 


The Anomaly B progran proviied a different oerspective 
from which to observe ta2 data sets since it prasented two 
dissimilar plots on the screen at the same time. The Lower 
traces consisted of five Anomaly A styl2 three-limensional 
Opets, this tim> only 2ijgaty points (29f bins) wide, while 
tne upper five traces were the two-dimensional time vs. 
amplitude plots of bins salected from the traces oselow. it 
was felt that tnis different perspective wouli 2nhance the 
ppomaly identification 2ffort, that it sea! not was 


disappointing. 


Qf the eighteen data sats analyz2i (3, 7, #16, Toe ee 0s 
Poe esos? so, 40, 41, 51, 53, 51, 80, 125, and 139), 
Anomaly B proved useful in only two instances: (1) to 
provide a graphic presentation of the three signal types 


Heed 2h the Statistical Aialysis Section of this daper; and 
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pete lOyebBase Scan 932, Base Freq 9.7 MHz, Base Bin 100 


Figure 17 — ATTENUATED SIGNALS EXAMPLE 
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(2) to illustrate one vary interesting phenomenon which 
could not be seen by the Anomaly A pr2sentation. AS shown 
in Fig. 19 this phenomenon consists of an undulation in the 
amplitude of the five upper traces. This "lazy sinusoidal" 
Variation appears to affact some signals more often than 
others. The term "lazy sinusoidal" should be usei with care 
Since it is not meant to imply that there is any strict 
periodicity to be found in the undulation or that the traces 
look exactly like sinusoldis. However, gany of the signals 
did avidence the more smooth (and, hence, sinusoidal) 


Varlation seen in trace on2 in the figure. 


Although only relativ2aly smooth waveforms are depicted 
in the figure, this effest was identified over a large range 
@2 Signal types including those which kayed off and 0on; 
however, it was much easier to see this phenom2non on the 
smooth signals. This variation was seen in every data set 
analyzed aithough it occured more frequently ani noticeadly 


in some than in others. 


The cause of tais undulation could be: (1) variations in 
transmitter output power levels; (2) recelver SYSt2n 
processing; or (3) ionosphera-related phenonena including 


Paraday Rotation, Sporadic-e& and fading. 
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FREQUENCY 


Base Scan 588, Base Bin 120 


Peewee 9o- AMPLITUDE UNDULATION EXAMPLE 
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Mee Doc USSLON 


In conjunction with th2 graphical analysis aoproach to 
identify Signal anomalies attriputable to. ionospheric 
properties, an attempt was made to glean additional 
information from the BRIGHAM data by means of statistical 
algorithms; however, it should be noted that no attempt was 
made to identify anomalies by this approach. @ather, tne 
programs MIN/MAX 1 and MIN/MAX 2 ware developed ia an effort 
to determine how oft2n and by what anount (aaplitude in 
decibels) both signal wavefronts (scans) ani individual 
Signals (bins or sets d€ bins) change. It was hoped that 
this new data might prov2 useful in the jievelopem2nt of an 
elcoOrithm tO adjust the myctciad sensitivity controls found on 
the swept-tuned receiving equipment us2i to acsguire and 
BecOrd th2 BRIGHAM data. Although this effort aas nothing 
directly to do with th2 signal analysis effort reported 
Previously, 1t did take considerable time ani effort tos 
Write the programs and cr2ate the desired data base. In 
addition, this evolutioa proved extremely h2loful in the 
overall effort 9f understanding signal characteristics, the 


primary goal of this stuly. 


For the purposes of tais inguiry, three different types 
9f signals (one that stayed "on" through most of the signal 
duration, one that stayei (Lee aerc olga hero st.OL Ehe. signal 
duration and oae that apo2ared to key on and off throughout 


most of the signal duratioa) were chosen fron data sets 
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which were collected at various tines of th2 day. The 
quotation narks are used (here and throughout th2 rest of 
the paper) to imply a g2neral trend noted over the life of 
the signal not some unvarying Signal characteristic. Figure 
20 is an Anomaly B graph plot presenting classic 2xamples of 
the three Signal types. [fhe "ON" signal (exhibiting that 
"lazy Sinusoidal" amplitude variation mention2d in thea 
graphical analysis sectioi) is trace 1, the "ON/OFF" signal 
is trace 2 and the "OFF" signal is trace 3. Traces 4 and 5 
are merely copies of trace 3 to fill out the plot. A table 
giving all pertinent jiata on the twenty-foac signals 
analyzed may be found in Pig. 21. Dist hes: Luc ees = t 
frequency refers to th2 base frequeacy of the sat when the 
recording was initiatel while TOI means the time Of 
maceCCeption for the data set. An attempt was made to 
analyze two signals fron each oof four different time 
categories (interesting from the ionosphere's standpoint) 
,those being sunrise, nod1, Sunset and aidnight. With one 


exception, set 53 - 2109LOCAL, this effort was realized. 


Be. STATISTICAL ANALYSIS PROGRAMS 


The first generation algorithm isveloped, called MIN/MAX 
1, provided for a scan-by-scan (down entire freguency 
bandwidths - cross sections of many signals - at fixed 
times) analysis of a variable number of scans witiin a given 
set. MIN/MAX 1 identifiei each of the peaks aad valleys 
along the scans, chose the highest peak and Lowest valley 
and computed the average peak and average valley heights. 
The quantization levels given oa the data cards w2cre changed 
to decibels and the amount of attenuation appli2i at the 
time of acquisition was added back in to giv2= received 
Signal levels accurately. MIN/4AX 1 proved to b2 of only 


Small value since signal wavefronts vice signals themselves 
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FREQUENCY 


Set 7, Base Scan 200, Base Bin 40 


Figure 20 - THREE MIN/MAX SIGNAL TYPES 
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were analyzed. A program listing for MIN/MAX 1 may be found 
at the end of this paper. 


The next generation effort, tabbed MIN/MAX 2, Focused on 
specific signals (the three varieties delineated in the 
discussion section abov2) and adii2ed a number of mors 
valuable Statistical features. Lik? its predecessor, 
MIN/MAX 2 identified each of the peaks and vailsys, chose 
the highest peak and lowest valley and computed the average 


peak and average valley heights. 


Its next feature was the deternination of two different 
types of signal changes: (1) the comparison of sample point 
SUplLcCudes “Wlth a Slidluagj ceferace point t> determine 
long-term changes in signal levels; and (2) th2 comparison 
of adjacent sample point amplitudes to determine short-tera 
changes. The former was initiated with the bin 41 point as 
the reference point and all subSeguant points conpared with 
it until an amplitude was found which differed from the 
reference by a previousiy agreei upon number 9f decibels 
(for this test, 4, 6, 12 and 18 dB changes were used in 
Succession). That point then became the referenc2 point and 
the process continued for the entire 2.4 minute length of 
the signal. The number of such Long-term changes and the 
elapsed time between each change were recorded. This data 
weeurd then be used t9 attain the goal of this 2ffort, the 
gauging of how often and dy what amount the varidus 9 signal 


types changed. 


A further refinement 39£ the procedure was th2 creation 
of amplitude pAistogranrs of the 2lapsed tines between the 
Changes to illustrate the effect that increasing the level 
Seeeaeecipel Eluctuation which constituted a change had on 
both the frequency of changes and the distribution of 
elapsed times. Since the receiver sampling rate was 25 Hz, 


the smallest elapsed time possible between K 4B fluctuations 
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was 0.04 seconds. A bcief analysis of the First set of 
elapsed times led to tha division of the time line 
(abscissa) into twenty-one time bins with th2 inclusive 
boundaries given in Fig. 22 and labeled N1 through N21. As 
a further refinement, oins N11, N2 and the first segment of 
N3 were divided into individual 0.04 second bins labeied 
foie stnrough Ni-/7;e02-1 through N2-3 and N3-1 respectively. 
The two histograms gave aa excellent perspective From which 
to watch the changes in signal characteristic as the amount 
of decibel fluctuation tested was stepped from 4 t9 6 to 12 
aid, tinally, to 18 dB. 


iiewE ina OULD OL MINAMAX 2°°was, an addition to the 
peak and valley statistics aentioned above, a listing of the 
number of long-term changes that took place in 2ar2 of the 
thirty-seven tim2 bins of interest, the percentage or the 
number of changes under each time Din as compared with the 
total number of chang2s (this allowed for a more valid 
comparison of the over-all change picture when an increase 
ln decibel deviation ceqguired resulted in fever changes 
tabulated in the bins), the total number of long-tearna 
Changes (callei NUM in tae output) and the total number of 
Changes in adjacent sampl2 points (the short-tera changjes 
mentioned above). In addition, the mean ani standard 
deviation for the elapsed times between long-term changes 
fere computed. aA listing for MIN/MAX 2 may b2 found at the 


end of this paper. 


Ce. OBSERVATIONS AND CONCLUSIONS 


There were a number of surprises which result2d from the 
Martial computer output From MIN/MAX 2. Prost bene. sign 
numoer of both long- andi short-term changes ideitified was 


not expected. Out of the possible 3597 shanges (there were 
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at most that many scans, 2%f sample points, available for 
each signal of interest) the number of changes noted in some 
cases reached as high as 20900 foc the low2r decibel 
fluctuation levels. Second, the elapsed time histograms 
tabulated showed a marked tendency to be heavily damped 
exponentials in nature with the 0.04 second bin (N1-1) 
having the highest numb2r of occurences throuyhout. And 
third, there was a high degree of overlap between tha 
Statistical Characteristics generated for the three signal 
types chosen indicating that ther2 was no absolute way to 
Pano TOniZe aAesignateas "ON" MON/Orr"® or “OFF werely by 
mooking at iitS Statistics. What was expect2i (and did 
occur) was that, aS to2 number of decibels used for 
comparison increased, the exponential distributioas "flowed" 
to the right indicating jceater elapsed times between larger 


decibel fluctuations. 


Figures 23 - 29 represent the bulk of the meaningful 
data garnered from the MIN/MAX 2 application to the 
mwenty-fLour signals of interest. In speaking of the 
Saabacteristics which differentiate the three irbitrarily 
Siosen dnd defined signal types, Figs. 23 ~ 28 provide those 
Peaelsottes | Which Cam b= Of valuS in coming t> a Logical 
GOnclusion concerning their differences. As would be 
expected from "ON/OFF" signals, they showed 4a greater 
difference between both average peak and valley values and 
between "maxmax" and "ninnoin"™ values than did the other two, 
lmplying more Signal leval fluctuation. The "OFF" signals 
consistently showed th2 smallest fluctuations as would be 
sepected from a series of either highly decayed signals of 
high but flustuatioajy aabient n91iSe.’ Basei on these 
meatistics alone it woald be nearly impossible ae 
distinguish with a high jiegree of certainty betw2en an "ON" 
and an "ON/OFF" signal jyiveno only the parameters in those 
figures. Figure 29 provides both a summary of the dat3 


concerning the number of shanges tabulated for 2ach decibel 
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comparison level and a general means of distinguishing 
between the three typ2s of signals. Althouga there are 
definite differences between the average figures for the 
various dB levels and some differences in th2 ranges of 
values, the overlap encouatered in the ranges makes signal 
differentiation by this means risky. In saort, what 
distinguished the "ON" fron the "ON/OFF" and "OFF" sagnals 
Was thoroughly subjective and Only loosely corrdborated py 


the statistics generated. 


Rddigvuenal informacion Can be retreived from Figs. 26 - 
28 which provide the nunber of changes (both long~- and 
Short-term) found when th2 decibel differential wis "walk2d" 
meome 4 toeiS GB and the oereentage of long-t2tm changes 
which fell into the most meaningful of th2 tairty-seven 
histogram bins (for elapsed tines between LONG] Som 
changes). AS 1s clearly shown, the number of caanges, NUM 
and ADJ, decreased as the dB level increased while the 
percentage figures in coiumns Ni and N1-1 2avidenced a 
movement away fron the ficst elapsed time bin. BOthH Of 
these characteristics imply that tae average time between 
Signal level changes at the various dB Settings iacreases as 
the number of 1B increas2. To be specific, by taking the 
data from Fig. 30 it could be calculated that 22 average 
time elapsed between 4 and 6 dB rluctuatuions foc the three 


Signal types were: 


Signal Type 4 dB 6 dB 
ON Sido oce Oe Sot 
ON/OFF 03.3 Ssocuwer 132° Sec 
OFF ~110 sec moe 9! Sec 
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Hirssdata GOnNELEMS taatethe "ON/OFF" signals fluctuated 
more rapidly than the other two categories while the "ON" 


Signals proved to be more stable than the "OFF". 


These results imply that, 1f£f the agency that operates 
the swept-tuned receiving system used to acquire and record 
the BRIGHAM data desices to acyuire FSK-ty9e2 signals 
Peoughly, the "IN" variety), 1t could gat away with changing 
the senSitivity adjustn2nts much Less frequently than if 
Morse~type signals (roujaly, the "ON/IrF" variety) were 
desired. Mito pe DLS tie mE EOuUgGH data provided la Fig. 23, 
Soula provide a good data base fron which a dial sontrolling 
algorithm could be adjusted to meet the needs of the signal 


types of interest at any tine. 


AS was mentioned above, MIN/MAX 2 provided for the 
calculation of the mean and standard deviation of the 
elapsed times tabulated foc eaca of forty-eight signals (the 
twenty~fourc signals of interest at both 4% and 6 aD 
comparison dJlev21s). Altasdugh the asan values wece accurate 
the standard deviation values were adt since the formula 
Moed in their calculation is valid for only normal (or 
SausSian) distributions. Since tne elapsed tine histograas 
were apparently exponential in nature, anoth2r form of 
variance determination was reguired. This effort called for 
Peepcucve-Tittingd technigus to identiry the histograms as 
being of a particularc distribution type, us2 of thea 
chi-squared test for goolaaess-of-fit and manipulation of the 
chosen distribution's meaa and variance formulas to derive 
the information desirei. Attempts were mad2 t> Fit three 
distribution funstions ts the histograms (Samma, 2?31SsSOn and 
Geometric) with the Sanna function proving t) be the most 
successful. The Gamma dastribucron Le come With 
parameters Ces, DOO <>, 0<A <= is defined as Follows: 


98 





f(x) = 


0, x<0 


: a 
Where, Mean -a%2, Vaciaace =AP . 


An algorithm was written which ceai in th2 observed 
histogram frejguencies, calculated the expested Gamma 
freguencies (using an International Mathematical Subroutine 
Library subroutine called MDGAM) and calsulated the 
chi-squared values (essentially, the sum of the least squared 
differences between the observed and expected frequencies) 
for use in the goodness-of-fit test. The chi-squared values 
was then compared with the appropriate value found in the 
chi-sguared distribution table (appropriate as to degrees of 
freedom and significance level); if it f21l below that table 
entry, the fit was considered good, if not the fit was bad. 
Romer scussed im Refs. 7 ani 8, a jg30d rile of thumb 15 that 
the smaller the shi-sguar2i value for a given nunber of 
degrees or freedom, the petter the fit. As a word of 
caution, even an 2xcalleat fit (extremely small chi-squared 
value) does not guarantee that the distribution is really 
the one being tested, it nerely states that the data do not 
present sufficient evidence to contradict the hypothesis 


that the histograms possess Gamma distributions. 


For the curve-fitting at hand, knowledge of tie mean and 
an iterative variation in xX  £4peraitted an 1terative 
Seeettation of 6 thus providing ail of the necessary 
lnftormation for tae algocithm discussed above. [a2 process 
was refined until the smallest chi-squared value could be 
identified. It should be noted that only the N1 through N15 


time bins were used since very few alapsed tines aver fell 


S)e. 





beyond five seconds and that, as the individual histogran 
compoSition dictated, either seven or three bins were 
actually used in the algorithm (following a suggestion that 
bins containing five or fewer counts be combined to enhance 
the accuracy of the test). The bin limits impliei that the 


degrees of freedom used were five or one respectively. 


Or the forty-2ight histograas t3> be fitted, Six had 
counts in too few bins to be useful (would have inplied ze2ro 
degrees of freedom - not found in the taole), 2ignteen were 
pad fits and the remaining twenty-four provided good fits to 
the Gamma distribution. Figure 30 gives the results of the 
curve-fitting effort. Knowledge of “% and th2 n2an implies 
knowledge of 4 which, i2 turn, implies (by definition) 
Knowledge of the vari3ace. A combination of ta2 mean and 
Variance values should be useful in the determination of how 
often the sensitivity dials on the rec3iving system shouldi 


be adjusted for the desired signals t9 be acguirei. 


Figure 31 provides a breakdown ofr the & values ys. 
Signal type and amplitide change Variable LOE the 
twenty-four histograms which "fit" the Gamma distribution. 
femewas frelit that this @*Would perait the JjJroduping or 
histograms with respect to some combination of the three 
categories above. The 31x ranges shown in the figure imply 
that any grouping by cg alone would be tenudus indeed as 
they overlap guite heavily. Thea only discernaple trend is 
the apparent lowering of the average M&M figice as the 
amplitude change variable increases from 4 to 5 dB. This 
would seem to contradict the defining equations for the 
Gamma distribution which inply that a lower & implies a 
more "damped-exponential" form for the distribution (when, 
in fact, the mor2 damped case is found at the 4% 43 change 
variaple). However, the shape of the distributioa is also a 
function of @ which, as iefined in Yean =%O , tends tod 


merect the distribution shape in a nore radical fashion 
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than does QA. In fact, as the change variable increases, 
the mean increas2s (se2 Fig. 30) ani & seems t) decrease 
(see Fig. 31) impiying 4 very pronounced ilacreasein Z 
wen, an the dsFining e27uations, is found t)5 be most 
prominent in the expoa2aeatial term; but, this implies a 
distribution shape less heavily dampled in the lower (4 4B) 
change variabl2 level. Paws Elena ot decreasing values, 
then, waS to be expect2i and appears to be the only 


meaningful product of the analysis of the Fig. 31 data. 


AS an aSide on the sadject of signal continuity during 
the brief recoried tim2 span of th2 waveforms (around 2.4 
minutes), it should be noted that data set number 7 bins 81 
ClON/OFF") and a (COL E*) exhibited VOcry unusual 
characteristics when compared with the other signals in 
theic classes. VogsvOmmm One lads c yon 8) Stil Pehadeoyer 
1000 changes at 13 1B difference while the other "ON/OFF" 
Signals had fallen off td much smaller numbers. In fact, 
that bin had two changes at 54% dB when ever-increasing 
numbers were applied. Th2 "OFF" signal, bin 91, 2videnced a 
Similar degree of consistency as it had many f2awer 4 JB 
Shanges than the other "OFF" signals while falliag off much 
more qguickly as the number of decibeis increased (see Fig. 
23) . In fact, this sijnal had only a few more than 1700 2 
dB changes implying exceptional amplitude continuity when 
compared with amy OF the twenty-three oth2r signals 
mialyzZed. TO giVe 2 NOES grapnic cepresentatio. of this 
phenomenon, Figs. 32 an41 33 compare the long-term signal 
Changes round in data sets 3 and 7 with the lLatt2r proving 
eo De the more consistent of the tv5, although Sven set 3 
mes 2 JOOd Measure Of unifornity throughout its 3600 scans. 
The base scans for each of the five scan sets (whicn show 
the beginning, nididle and 2nd of each data set - time-wise) 


Mee provided at the right of the pictures. 


62 





TIME 





S06 


Es \ t NaN \ . S / \ \ - 
ee y /, ‘ J \y wo iN Y ees ENT a i 
fas ‘) \ oa, oF x | Vea le \ i\ ~/ a i, 
O : ae \ ae ae = | AEZ \ ! E / ~Y/ ‘ y, u 1 6 


/ 
\ v fe ee W/ 
\ a e ea ees ee f \ awe \ vs ‘ vx / 
) af ai iV. NA INAH, \ / ce a Be: ae oe 
i Va } \ / a ea 
in RAR LRS 
‘ ey 7 \n \ ay. wa a ANS, ry Se 1816 
/~ a 4 j 
“\ aaa ONIN at 
“ia ; wy las _\ - ' \ yas 
Degen TAWA aA eo fe 
ax oy) = o- : Bt zo Ue x 7 San \\ a are 
ae aera AMAT Son AS 
YARN MA VE VAAN ONS 2226 


ee a aX. | \ Dp 8 ae ef 
Nn ~ SA Wi PS AY . \ va cc iy 

oe y; As \ Ne ON pel 
wii Awl u Nav SPARE ) re 


/\> [\- / IN A 4 
MIRAE SA YL Ce POF 





FREQUENCY 


Paige Ted aovoeMHe Base Bin 100 


ieee 2) - LONG-TERM SIGNAL CONTINUITY = SET 3 
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LONG-TERM SIGNAL CONTINUITY - SET 7 
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V. COMMENTS AND RECOMMENDATIONS 


Oeetie two, Signal anomaly iientification tEechnigues, 
Anomaly <A provei to be the most fcuitful. Naccow and wide 
PimceerOveeatScCONCInUlties, <fLading, frequency Shifting of 
entire wavefronts, frequency shifting of individual signals 
and several miscellaneous varieties were seen. [wo of the 
above (discontinuities 2121 shifting) «ere reported on in 
pete | Wath the rest being unigue t> this study. Several 
examples of pseudo-anomalies were includjed to gerform the 
dual function of warning against the false assimption of 
anomalous perfornance sinply because the plot seemed to 
1mply it and to provide pictures with th2 narrativ? in such 
cases. Analysis with Anomaly 8B proved to be celatively 
unproductive as only th2 iaterasting and freguentiy seen 
phenomenon of signal amplitude undulation (reported in Ref. 
1 and expanded upon herei= oveing noted. The diverse and 
amazingly compl2x natace 9f the iLonospher2 and its 
interaction with its e2ivironment precludes tn2 precise 
assignment of a cause td every effect (1.2., anomaly). 
Suffice it to say that, except for thos2 anomalies which can 
be explained away by eguipnent malfunction, nothiag has bzen 
meuma If this Ssf££LIrt whica will simplify the modal of the 
HonoSphere aS an ever-Shanging anelium in whici literally 


anything can happen to a £ad19 signal which travecses it. 


The statistical analysis work performed, from creation 
Meeene MIN/MAX algorithns ts the tedious surve-fitting 
effort, resulted in 4a wealth of predominantly raw data 
Which, with the mean ani vaciance values, shouli provide a 
Foundation upon Sook aA ceceliver systen sensitivity 


Meaet-cOntrOlling algoritnn could be refined. [t was nota 


in the curve-fitting section that it is perhaps the weakest 
Ihemietne tae Statistical study Since the chi-sjiared test 
guarantees neither that the elapsed time nistograms are 
Gamma distributions nor that thera tight 2xist other 
fester puttons whlch Wolzld “E£1t" the histogram data as well 
as the Gamma function iid. However, a Statuse ve ql isy 
reasonable argument can be nade that the histograms (at 
least half of them) ars probably Ganma in nature and, hence, 


the variance values calcialated can b2 accepted as accurate. 


Should additional work be done with the BRISHAM data, 


the following reamarks/rceconnendations could prove useful: 


(1) all of th2 compleat2 data sets in hand have been 
thoroughly analyzed witn program Anomaly aA - additional 


effort along this line nigot D2 a waste of times and energy; 


(2) althougno a2early twenty of the data sets in hand have 
not been analyzed with Ddrogram Anomaly 8, the disappointing 
results thusfar might act aS 4a warcning as to axpacsted futures 


Bewalds from adiitional effort; 


(oy) 2f AnoOMaly A and/St Anomaly 8 ars to be ased again, 
the BAKSCN featurs found ion tha METASYMBOL subroutine (it 
provides for tn2 backiag over of a Specifiei number of 
files, or data sets) should be nade to work aS a time/labor 
Saving device - the problem 1S in the magnatic tape unit, 


mot the progran; 


(Qi rae — EFscOMmMendiation (nade initially in 222. 1) f€9e 
mae Creation of a computer algorithm t95 seatcn out the 
anomalies in an automatic and, hence, more rapid aanner than 
the visual approach 1s jgaacdedly endorsed herein - although 
the speed at wnich the data sets could be analyz2i would be2 
improved by at least an order of anagnitude, t12 extreme 


Ponplexity of trying t> tell the sconputer what to look foc 


66 





as well as the loss of that very intangible but valuable 
human discretionary power sould combine to made this effort 


a quagmire; 


(5) should additional work be don2a utilizing the MIN/MAX 
routines, Anomaly B shouli pe ased first to eahance the 
peeclracy in selecting tos different Signal types to pe 
analyzed ~- that this was aot done in the present study aay 


have resulted in the rath2ar Large ranges found in the data; 


(Vereen naCcCuULAGCY Of Bly LUtarpe Subve-f£itting work coulda 
be enhanced by the inslusion of the Weibull and Lognormal 
SeStELoUtions (96 day 29f the Jlevel-srossing distributions) 
to the list of possibdl2 curves td whith the selapsed tiase 
histograms are matched b2acaise sets of such 212ata can be 
fitted to more than one Tistribution with equal ascuracy and 
validity - this was not jiodne here because th2 suggescead 
PeStELOUtCIONS were complex e2n5ugn and the tins remaining 
Short enough that a good 2ffort couli not have oeen put 


morth. 
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GLOSSARY 


Poorer ite tec emonIZontal layee {OF patcii O15 oT Pilouars 
ay2 


V2 


electron density enbedded in the regularc E | 


BESipbe rng In Signal deflection and/or absorption. 


Pane alaye AsyiMeatry: ~LlassS attenuation at nigat than day 
gieweordirrerence in Ehe Lonizingj effect of the sun 
resulting 12 discoatinuities iaring the sSuncise/sunset 
periods. 

Geomagnetic Non-Resiorocity: different attenuation 


effects on signals travelling wast to east rather than 


east to west due to the earth's nagaatic field lines. 


Pepese eet ADSOU DELON (PEA); absorotion of FradiS signals 
Pascoe nguprge= —S0lar Eegions = dus te a comDlLnation of 
solar flare activity and the shaose of the earth's 


Magnetic field lines. 


RUEOudImDL So Maye — bDrigatly scoloced aodrthecn latitudis 
phenomena disruptiv2 t9 coOmmdnications - also as the 
moSUineOrmt  COMbInation of solar flare activity and the 


earth's magn2tic field lines. 


Fataday Rotation: any Lineaarly odlarized wave travelling 
iene beect ion of 4 Magnetic ELilei results ln Lets tw 
Circilaaly polacizsei sanponents Erayvealling at different 
velocities and thas the plane of polarization will 


hOtate WhLth distanc:. 
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Peeevilcelgadssar a Fadio wavs comes int contact with a region 
of inhomog2neous refractive indices the wave could be 
trapped between twd layers ani guided, as ina leaky 


wave guide, away fron its intended destination. 


8. Fading (Dellinger Effecty: sudden ionospheric 
disturbance (S. I. 9.) which proiuce2s a conplete radio 
"fade out" jdastingjy from a few minutes to an hour or 


more - caused by solar flare astivity. 
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PEOLTeROULTINE. “aa.ec 


Hdefine  XINT 320 
tdefine YINT 8000 
Hdefine DELY 3200 


A#define DELX 16 

H#define AMPY 8000 

A#define NLINE 1 

tdefine PI 3.141596525 
Hdefine NPT 150 


int zero 0; 
int one 13 


Struct { 

int not; 

int line; 

int Idx; 

int Idy; 

int cotrs 

int nspt? 

int spotl: 

int spote? 

int spot3; 

int hremineg sec? 

int nft, 

int laa; 

int ser; 

int fd(3); 

int Ip; 

int scl; 

int mordayryr: 
) head, th; 


int 4buf (300); 


Struct data ( 
int x3 
int yds 

)3 


Struct data *d; 
int dbuf (4000) 3 
int np? 

IiNt nip 8B; 


int idevre,odev.esclxsscly-hiasx,tdev; 


main(arac,argv) 
Int **arav; 
{ 
iNt jejons 
char *¢s;7 
if((tdev = open("/dev/spp",1)) < O)4 
printf("“cannot oven spon 4"); 
exit(); 
} 
if((nmdev = open("/dev/rvn”",1))<0){ 
printf ("cannot open ryn 4"); 


L00 





} 
ne Gl 


} 


exit); 


idev = open("/dev/rmt6",0))<0) { 
printf("cannot open rmt6 4"); 
exit); 


if(araqe > 1){ 


} 
whil 


stty 
} 


#define 
#define 
k#define 
4#define 


Struct ji 
int 
char 
char 
char 
int 
int 
are 
int 

} itabt7 


char nbf 


cs = arav(l}; 
n = 0; 


while ( *es >= °O* RR tes <= '9*) 


m=n tr 10 + *estt - *0°? 
n= 113 
for ( 420% i<ns itt) 
inn(idev, dbuf,800); 


orintf ("number of records skipped, 


e((CinplCidevrhbuf,809))>0) ( 

h = hbuf? 

j = 0; 

for(i 0; i< HK => Lines +4) 


jo =+ e*(h->not);? 
} 


scilx = 0303 scly = 033 
biasx=50; 

Np = h=>nnt * h->line; 
scale(); 

plot); 
cvers(pdev,020); 
(pdev,r&one)? 


NBYT 264 
NBLK 0 

NSL 1250 
DRAW 1 


pt { 
Ye 
*xoh; 
yinc; 
xdir; 
yf; 
cxnes 
EV VM ¢ 
thlink; 
ANB V1) 5 *10,7%* tas 


NBYTI; 


EO ae 


n = inplidev,&dbuf (jj), 800); 


Zd#",n)5 





int *dp[2000); 


plot () 
{ 


struct tpt *s; 
int irjries 


sort (); 
for(i = 0% + < NBYT; 449) //clear olot buffer 
po{i) = 0; 
for(i = OF 31 < NBLK; +4) 
write(pdev,obs/NBYT); //move to top of plot area 
ip = 0; fa = itab; //plot scan line 
j = 07 
for(i = NSts £ > 1s ten) { 
while(j < no && *tdol(j) =+ j)f 
ie = sio(dplj}); //set up plot point for interpolation 
if(je == 1) 
return; 
J++; 
} 
nib(id; //set up Plot buffer 


write(ndevr,ob,NAYT)://plot line 


int jsort? 


sort() 
{ 
register iskr»t? 


1-0; 
for(d = dhufs d < &dbuf [nore) > dt) 
apotit?)] “= &l(qg > yd); 


k=np; 
while ( k =>> { jf 
Jsort++; 
while ( jsort ){ 
jsort = 0; 
for € £207 §< (npek); it) 
me CMe delitns, «tdoli+k) ) 4 
t=dnofi)}; 
dol(ilsdpolitk); 
doli+k)=t; 
jsortt++; 
} 
} 
} 


2 


siop(do}) 
int *dpi; 


{ 
int levi kis Via k Pr VPs WNC Vy tVrdirxrexr*S- lox, inc; 
int tries 
s=dpi; 
lex = I? 
yi = sedptrer;s xi = kdoiz 
if(++s < &Bl(dbuf (etnpl] )) 
xP = *gtt; 
yr =' ks; 
if€yr & DRAW) { 
if( Cine = yi - yr) >= 0) { 
if( Cincy = inc) == 0) 
dirx = nip; 
else { 
for(i = 03 (Cincy = Cine/(nip >> 4+))) =s 
dirx = (1 << 4);3 
} 
fy = yr; 
cx = 0200; 
lox = 0; 
je = stacklyvicrincyrincyrfyrxiredirx,cx)? 
Gass 1) | 
return(1)3 
} 
} 
} 
if(eedpti >= dhuf) { 
if(yi &8 DRAW) { 
yl = t*dpies; 
xl = kdpi; 
if( Cine = vi - yl) >= O){ 
if(Ciney = inc) == 0) 
dirx = nip; 
else 
fortis = 07 (CCUimcey = Limne/(nip >> i))). == 
Girx = =(1 << 4)> 
} 
fy = yl; 
cx = lex; 
je = stacklyirincyrincyrfyrxirls,dirxr,cx); 
veCie.== [) 
return(1); 
} 
} 
} 
return(0)3 
} 


Stack(a-sb,cr,dd,e, f) 
int arhbrec,dd,e,f; 


OS 


0) 


0) 


e 
s 


e 
v 


1t4) 3 


Tet 





{ 


//setuno olotting 
//delete point 


//Yine break 


int *3;7 
int 2 
Struct tot *27 
ja-> y = a; 
ia => yine = oO; 
ia => yf = ce 
if(dd >= NBYT && dd < 0) 
printf("bad scale x«=%d4#",dd)i 
exit(); 
} 
13a °-> xpb = dd + pb, 
Poe => ex Gt ft -=.<e> 
ia => cxop = fz 
if(iop == 0) ¢ 
ip=itab; 
{op->flink=to-=>blink=0; 
iatt+; 
iaz>blink=ip; 
ja->flink=0; 
} 
else{ 
Sz=ia=>blink: 
s*>flinks=ia: 
V#Cra->flink == 0)4 
S$ = {ar 
if(ttia >= Litab(7*NBYT) ){ 
printf("itah overflow 4"); 
return(1); 
) 
iar->blink=s;? 
ja=>flink=0; 
) 
eise{ 
s = ia -> flink: 
$s -> blink = ia? 
ia -> flink = 0; 
ia = S?> 
) 
) 
) 
nib(sl) 
int sl; 
( 
int ir jrne 
int *S7 § = 103 
while(s) ( 
if(s -> yf <0) 
if((szfree(s)) == 0) 
return; 
Sa ee OO. = Sa => CXD sr 
if(sl == $s -> v)f 
ye= sC(mes G => edir)y > 0 7 
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n 


= ss 


buffer 





if( n < 90 ) //i\eft 
fort) = 03 4° <- Va jee) 
14:05 “*>¢xp. == 0) 4 
#g => xpb =¢ 1; 
$s *> cxp = 1; 
} 
ts => xob =: (5s |-> cxp =<< 1); 
} 
else 
for(j = 0% 3 < f2 jee) 
ager>xnb =) (s->cxp =>> 1)3 


} 
S*>y =- s*>yINc? 
} 
if(s) <= s->yf) /fend of point 
s->yf = ol; 
§ = se>flink? 
} 
} 
free(s) 
int *3$7 
{ 
int *t; 
int i? 


struct int *2; 
5 => xypb = ‘ 
if(s-r>blink == 0){ 

Di c= ~s=>f link: 
{p2=>blink = 0; 
t=1p7 
} 
else{ 
t = serblink? 
t~e>flink = g2r>flink: 
t = s => flink; 
if(e ss 0) ( 
s => flink = fta3 
ja = 83 
return(t)3 
} 
t => blink = s 2> blink; 
} 
s*>flink = ja? 
s->blink = fa->blink; 
ja = 83 
return(t); 


scale() 

{ 
Struct data *s; 
iINnt Omeiv fr Ce MINesMINV: 
int dxl20) ,dy(e9); 


OS 





} 
inp ( 


{ 


} 


char 


conv 


$ = dbuf;> minx = miny = O77717; 
d=s, 
for(i = 07 31 < nope 144) 4 
minx = ((¢ = § => x) < minx? ¢ 3 
miny = ((€c¢ = § *-> yd) < miny? ¢c 3 
S++; 
} 
if(minx > 0) 
minx = 03 
if(miny > 0) 
miny = 0; 
s=dbuf; 
for(i 2 03 4 < nps i++) { 
Om = § => yd & 13 
Ss => x ==> Minn; 
$ => yd =" miny?: 
62-55% =/ “Sclx; 
S$ 7>x 2+ btasx; 
af (se>x >= NRYT) { 
printf ("overflow po #4"); 
exit(); 
} 
$ => yd =/ scly; 
if(s=->yd > NSL) 
s->yd = NSL-13 
$s => yo =% 0177776; 
$ => yd =; dm; 
Stt,; 
} 
idf,buf,nbyte) 
int idf,*buf,s,nbvte; 
int frtsysnecs 
struct ( 
char clrcerc$,c¢ 4 
» cff{1200), *s; 
$s = cf; 
if((n = readlidf,cfsnbyte)) > O){ 
for(i = 07 4 < nbyte/43; i++) { 
Go =. 5, hee <<-.e; 
oso 1, 
t=; § => ¢c3 << 6; 
t =; $ => ¢4; 
S4#+; 
buf fil = t; 
) 
} 


tcbop,tloc,tllise2l -ch lio}; 


(val) 
int va): 
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minx); 
miny)?; 





int a? 
if€(a = val/10) 
conv(a); 
febort = yal -% 10 + '0"s 


conc(cl,ce.n) 
char *c23 
iNnt mexrel? 


(- 
int jom; 
fon CUrve= 0708 <n <}44 ut 
1f(i == 0) 
m= (tc! 2% 07700) >> 63 
else 
m= *¢€1 & O07 7; 
itim == O02) 
tee = 0 
if(m >= O01 &% m <= Ott) 
See S 9 | tm = ]7 / 
if(m >= O21 && m <= 031) 
koe = 'AT + m ~ O2T3 
if(m >= O41 && m <= OST) 
tee = 'J' + m -041;3 
if(m >= 062 && m <= O71) 
zee = 'S' + m = 0623 
cett}; 
} 
return(i); 
} 
skion(ent ) 
int cnt; 
{ 
nt 44 
for(i = OF | < ents 14+) 
write(odeverph,2); 
} 
clr) 
{ 
int i? 
fortil= OF) <1 322 744} 
Cris =" 3%; 
tP(i3st)] = *#*; 
} 


mov(to-rfrom,n) 


Oy 








int ne 


int i? 


forty = Om < ae 44") 
etot+ = rf romtet? 
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